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Abstract 

The present study was performed to investigate tine variation of pliytocliemical composition, antioxidant activity and High 
Performance Liquid Chromatography (HPLC) fingerprinting profiles of three Pyrola species. Thirteen samples (eight P. 
decorata, three P. calliantha and two P. renifolia) were collected from different regions In China. The tannin, hyperoside and 
quercetin contents of all samples were determined by reverse-phase HPLC and varied within the range 9.77-34.75, 0.34- 
2.16 and 0.062-0.147 mg/g dry weigh, respectively. Total flavonoid content was evaluated and varied within the range 
16.22-37.82 mg/g dry weight. Antioxidant activity was determined by DPPH assay, with IC50 ranging from 7.96 to 50.33 \ig/ 
ml, ABTS'^ and FRAP assay, within the range 612.66-1021.05 and 219.64-398.12 |xmol equiv. Trolox/g, respectively. These 
results revealed that there were significant variations in phytochemical profiles and antioxidant activity among all samples. 
Due to the higher phytochemical content and significant antioxidant activity, P. calliantha was selected as the most valuable 
species, and the P. calliantha sample from Left banner of Alxa even possessed the strongest antioxidant activity among all 
the thirteen samples. Futhermore, Emei IVlountain was proved to be the most suitable region for producing P. decorata. 
Moreover, in order to further evaluate the diversities and quality of Pyrola, HPLC fingerprint analysis coupled with 
hierarchical cluster and discrimination analyses were introduced to establish a simple, rapid and effective method for 
accurate identification, classification and quality assessment of Pyrola. Thirteen samples were divided into three groups 
consistent with their morphological classification. Two types of discriminant functions were generated and the ratio of 
discrimination was 1 00%. This method can identify different species of Pyrola and the same species from different regions of 
origin. Also, it can be used to compare and control the quality of Pyrola and other natural products prepared from them. 
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Introduction 

The evergreen herbs of the genus Pyrola are traditional Chinese 
medicinal plants. Pyrola is mainly distributed in the northern 
hemisphere in temperate and cold temperate regions around the 
world with about 30 species. In China, there are 27 Pyrola species, 
mainly distributed in the west and northeast of the area. [1]. Pyrola 
has very high nutritional and medicinal value: it has been widely 
used in functional foods and supplemental products for treatment 
of rheumatism, waist pain, knee pain and high blood pressure [2- 
4]. The water extracts of the plant were reported to inhibit the 
growth of many kinds of human pathogenic bacilli in vitro, which 
were used in refreshing foods [5] . Nowadays, because of its health 
benefits, Pyrola is consumed as a kind of tea called "Lu shou cha" 
for daily drinking in China [6]. 

Pyrola contains a range of bioactive components, such as 
flavonoids, phenols, quinones, terpeiies, amino acids, etc. [7,8]. 
Polyphenolic compounds are especially abundant in Pyrola, which 
may contribute to its high antioxidant activity. A previous study 
has reported that the ethyl acetate fraction from P. incamata leaves 
possessed remarkable antioxidant activities, which were almost 



comparable to the capacities of vitamin C and butylated 
hydroxytoluene [6] . Another study has also found that the radical 
scavenging activity of P. incamata from Tahe was very close to that 
of vitamin C [1]. Taken together, these previous studies indicate 
that Pyrola exhibits powerful antioxidant activities and has 
potential application in nutraceutical and functional food products 
for improving human health. 

The chemical compositions of Pyrola vary gready since it has 
various species and extensive distribution. A previous study from 
our laboratory demonstrated that ecological factors contributed 
significandy to the contents of active compounds in Pyrola samples 
[9]. Another study from our laboratory indicated that the metal 
element contents of P. decorata varied immensely from region to 
region [10]. In addition, Zhang et al. reported the great variation 
of active components in P. incamata from eight sites in northeast 
China [1]. Different chemical compositions may lead to significant 
differences in the effectiveness for health problems and the safety 
related to intake [1 1]. Therefore, an evaluation of the diversity of 
Pyrola from different sources would be desirable in order to ensure 
the quality of Pyrola and its derived edible products. 
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The objectives of this study were: (1) to analyze the variation in 
phytocliemical composition and antioxidant activity (DPPH, 
ABTS and FRAP assays) in three Pyrola species (P. decorata, P. 
calliantha and P. rmifolia) collected from different regions in China; 
(2) to obtain the HPLC fingerprints of Pyrola and establish a useful 
method for identification, classification and quality assessment of 
Pyrola samples using HPLC fingerprint analysis coupled with 
hierarchical cluster analysis (HCA) and discrimination analysis 
(DA); and (3) to provide meaningful information for the selection 
and application of Pyrola in both healthcare and the food industry. 

Materials and Methods 

Instrumentation and reagents 

HPLC analysis was carried out with an Agilent Series 1260 
liquid chromatograph, equipped with a quaternary gradient pump 
system and a variable-wavelength detector system, connected to a 
reverse-phase (RP) SB-C 18 column (5 |a,m, 4.6x250 mm, Agilent, 
USA). Data collection was performed using ChemStation software 
(Agilent, USA). 

Hyperoside and quercetin were purchased from Sinopharm 
Chemical Reagent Co. Ltd (Shanghai, China). Tannin was 
provided by the Tianjin Dongliqu Tianda Chemical Reagent 
Factory (Tianjin, China). Chromatographic-grade methanol and 
CH3CN were purchased from Tianjin Bodi Chemical Holding 
Co. Ltd (Tianjin, China). Other chemicals were of analytical 
reagent grade and were purchased from Tianjin Bodi Chemical 
Holding Co. Ltd (Tianjin, China). All aqueous solutions were 
dissolved in deionized water. Stock solutions of tannin, hyperoside 
and quercetin were prepared in methanol and were diluted to the 
desired concentration. All solutions were filtered through 0.22 (im 
nylon filters before use. 

Plant materials 

Whole herbs of thirteen samples (eight P. decorata, three P. 
calliantha and two P. renifolia) were collected from diflferent regions 
in China (Table 1). All materials were air-dried and powdered, 
and were stored at — 20°C and protected from light until further 
analysis. 



Ethics statement 

Specific permissions were not required for the described field 
sampling studies or for the collection of plant materials. For any 
locations/activities, no specific permissions were required. All 
locations where the plants were collected were not privately owned 
or protected in any way and the field studies did not involve 
endangered or protected species. 

Preparation of the extracts 

Each powdered sample was extracted twice with ten times of 
95% ethanol for 2 h in the dark, under nitrogen, at ambient 
temperature and pressure. All extracts were stored at — 20°C in 
the dark for further use. Extracts were diluted if necessary. AU 
extractions were performed in triplicate. 

RP-HPLC analysis 

Sample solutions were thawed, filtered through 0.22 H-m 
membrane filters, and then separated by RP-HPLC to obtain 
chromatograms. The amounts of three polyphenolic compounds 
(hyperoside, quercetin and tannin) were quantified using RP- 
HPLC at ambient temperature. The mobile phase consisted of 
H2O-CH3COOH (200:1) (solvent A) and CH,,CN (solvent B). 
The flow rate was 0.8 ml/ min. The gradient program was set as 
follows: 0-8 min, eluent B was kept at 18%; 8-9 min, eluent B was 
increased from 18% to 20%; 9-17 min, eluent B was kept at 20%; 
17-30 min, eluent B was increased from 20% to 34%; 30-35 min, 
eluent B was increased from 34% to 65%; 35-45 min, eluent B 
was kept at 65%. The injection volume was 20 |xl, and the detect 
wavelength was 370 nm. Analyses were performed in triplicate. 

Determination of total flavonoid content 

The total flavonoid content was determined by reactions of the 
AICI3 method with some modification [12]. Rutin (4-40 mg/1) was 
used for calibration of a standard curve ()'= 13.067x+0.0367; 
R = 0.9996) and results were expressed as milligrams rutin 
equivalent per gram dry weight (DW) of sample. AU samples 
were analyzed in triplicate. 



Table 1. Three Pyrola species (P. decorata, P. calliantha and P. renifolia) collected from different regions in China. 





Species 


Sample Number 


Origin 


P. decorata 


SI 


Dianbingchang, Taibai Mountain, SX 




S2 


Lujuanliang, Taibai Mountain, SX 




S3 


Shanyang, SX 




S4 


Ningshan, SX 




S11 


Dadian, Taibai Mountain, SX 




SI 2 


Fangyangsi, Taibai Mountain, SX 




S7 


Emei Mountain, SC 




S9 


Xiabaiyun, Taibai Mountain, SX 


P. calliantha 


S5 


Zhongshansi, Taibai Mountain, SX 




S6 


Left banner of Alxa, NMG 




S10 


Huzhu northern Mountain, QH 


P. renifolia 


S8 


Xingshan, HB 




S13 


Linjiang, JL 



SX: Shaanxi province, China; SC: Sichuan province, China; NMG: Nelmenggu province, China; QH: Qinghai province, China; HB: Hubei province, China; JL: Jilin province, 
China. 

doi:l 0.1 371 /journa!.pone.0096329.t001 



PLOS ONE I www.plosone.org 



2 



May 2014 | Volume 9 | Issue 5 | e96329 



Diversities and Quality Evaluation of Pyrola 



300 -, 



250 - 



< 200 • 
S 

E 150 . 



5 100 - 



50 - — 



0 -I 




1 — ' — I — ' — r 

15.0 22.5 30.0 
Rentention time (min) 



1 

45.0 



Figure 1. HPLC chromatograms of Pyrola samples. Three 
compounds were identified by their retention times (min): tannin 
(2.34, peak 1), hyperoside (13.26, peak 4), quercetin (32.28, peak 9). 
doi:10.1371/journal.pone.0096329.g001 

DPPH radical scavenging activity 

The DPPH radical scavenging activity was measured by a 
slightly modified method [13,14]. Amounts of 2 ml of the tested 
samples (10-90 |ig/ml) and the positive controls (quercetin and 
rutin, 1-50 |J,g/ml) were mixed with 2 ml DPPH solution 
(200 ixM); the final concentration of DPPH was 100 (xM. The 
mixture was shaken vigorously and allowed to stand in the dark for 
30 min. The absorbance at 517 nm was measured with a 
spectrophotometer. All the tests and the controls were repeated 
in triplicate. DPPH free radical scavenging activity was calculated 
using the following equation: 



SA(%) = [l-(A--4)/^olxlOO 

where A^: absorbance of ethanol (2 ml) and DPPH (2 ml); A;, 
absorbance of tested samples (2 ml) and DPPH (2 ml); Aj. 
absorbance of tested samples (2 ml) and ethanol (2 ml). IC50 
values were the effective concentrations at which DPPH radicals 
were scavenged by 50%, and were obtained by interpolation from 
linear regression analysis. 



+ 1 <!' 



ABTS*^ radical cation scavenging assay 

The method of decolourisation of free radical ABTS*^ was 
performed according to Re et al. with some modification [15]. The 
ABTS» was prepared by mixing an ABTS stock solution (7 mM 
in water) with 2.45 mM potassium persulfate. This mixture was 
allowed to stand for 12-16 h at room temperature in the dark until 
reaching a stable oxidative state. For each analysis, the ABTS'"*" 
solution was diluted with pH 7.4 phosphate buffered saline (PBS) 
solution to an initial absorbance of 0.700±0.021 at 734 nm. This 
solution was freshly prepared for each analysis. For the spectro- 
photometric assay, 100 |al samples (100 (ig/ml) was added to 
3.9 ml of ABTS* solution and the absorbance was determined at 
734 nm. Results were expressed in terms of micromoles trolox 
equivalent per g of dry weight (|J.mol eq. trolox/g). All 
determinations were carried out in triplicate. 



Ferric reducing power (FRAP) assay 

The method of FRAP assay used was a modified version of that 
reported by Benzie and Strain [16]. The method is based on the 
reduction of a colorless ferric complex 2, 4, 6-tripyridyl-s-triazine 
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Table 3. Content of tannin, hyperoside, quercetin and total flavonoids in Pyrola 


samples. 




Species 


Origin 


Content (mg/g DW) 










Tannin 


Hyperoside 


Quercetin Total flavonoids 


P. decorata 


SI 


14.26±0.06'" 


0.86 ±0.03"' 


0.080±0.006'" 21.42±0.47"' 




S2 


9.77±0.02' 


0.46 ±0.01" 


0.070±0.008'" 17.53±0.72^'' 




S3 


11. 84 ±0.05' 


0.35 ±0.01' 


0.064±0.014' 16.22±0.48'' 




S4 


14.90 ±0.06''= 


0.55±0.01'''' 


0.062±0.004' 20.74±0.45''' 




Sll 


9.77±0.01' 


0.82 ±0.01"' 


0.063 ±0.003' 20.90 ±0.44'' 




S12 


16.48 ±0.08'' 


0.96 ±0.03' 


0.076±0.003'" 22.84±0.48' 




S7 


12.62±0.0l'" 


0.34 ±0.01' 


0.067±0.001' 37.82±0.54'' 




S9 


10.60 ±0.07' 


0.57 ±0.01''' 


0.066±0.002' 22.80±0.33' 




Average 


12.53±2.50'" 


0.61 ±0.24''' 


0.069±0.006' 22.54±6.61' 


P. calliantha 


S5 


3 1.84 ±0.24' 


1.79 ±0.07'' 


0.1 31 ±0.006" 16.66±0.6l'' 




S6 


34.75 ±0.05= 


0.46 ±0.04' 


0.147±0.007" 19.30±0.34'' 




SIO 


20.34 ±0.1 2'''' 


1.93 ±0.03"'' 


0.084±0.002'" 18.31 ±0.5l's 




Average 


28.98±7.62'' 


1.39 ±0.81'" 


0.1 21 ±0.033" 18.09±1.33'S' 


P. renifolia 


S8 


15.47±0.05'" 


2.16±0.07'' 


0.074±0.012'" 26.77±0.52'' 




S13 


13.12±0.1l'" 


2.09±0.03" 


0.090 ±0.003'' 27.75 ±0.47'' 




Average 


14.30±1.66'" 


2.13±0.05" 


0.082±0.01l'' 27.26±0.69'' 


Each values represented in tables are 
Values with different letters (a, b, c, d, 
doi:10.1371/journal.pone.0096329.t003 


-neans ± SD (N = 3). 

e, f, g, h) within same column were significantly different (P<0.05}. 


omplex (Fe'^'''-tripyridyltriazine) to a blue-colored ferrous form in 40 niM HCl and 20 mM FeCls ■ 6H2O solution, with the ratio 
"e ^-tripyridyltriazine) by the action of electron-donating antiox- of 10:l:l(v/v). For each analysis, 400 |a-l of sample solutions 
dants. The FRAP reagent included 300 mM acetate buffer (3.1 g (500 ng/ml) was added to 3 ml of freshly prepared FRAP reagent. 
j2H3Na02 " 3H2O and 1.6 ml C2H4O2), 10 mM TPTZ solution The reaction mixture was incubated for 30 min at 37°C in a water 


Table 4. Antioxidant activities of Pyrola samples. 








Species 


Origin 


DPPH ICso (ng/ml) 


ABTS ([imol Trolox/g) 


FRAP (^mol Trolox/g) 


P. decorata 


SI 


44.56±0.01^ 


612.66±1.26* 


256.07±0.39' 




52 


38.23±0.0l' 


698.52±0.32' 


223.20±2.3l5 




S3 


27.74±0.02°' 


735.43 ±1.65'' 


31 1.45 ±2.04''' 




S4 


30.56 ±0.08'''' 


758.05 ±3.88'' 


306.45±3.94''' 




Sll 


42.85±0.1l's 


643.09 ±2.54'' 


223.85 ±1.85^ 




SI 2 


50.33 ±0.01'' 


615.75±0.99' 


219.64±4.663'' 




S7 


29.91 ±0.03'' 


798.23 ±7.33' 


309.1 2±1. 47''' 




S9 


32.69±0.14'''' 


712.50±0.23''' 


298.09±0.67' 




Average 


37.11 ±8.16' 


696.78±67.94' 


268.48±42.08' 


P. calliantha 


S5 


9.02 ±0.03'' 


958.66±1.45'' 


376.52 ±2.98"'' 




S6 


7.96 ±0.04'' 


1021.05±4.21' 


398.1 2±2.87" 




SIO 


11.99±0.06'' 


974.20 ±0.69'' 


355.04±0.6l'" 




Average 


9.66±2.09'' 


984.64±32.48'' 


376.56±21.54"'' 


P. renifolia 


S8 


26.56±0.03°' 


828.60 ±9.54' 


3 16.45 ±5.32'' 




S13 


23.03 ±0.02' 


809.32 ±1.27' 


335.41 ±1.47"' 




Average 


24.80 ±2.50' 


81 8.96± 13.63' 


325.93 ±13.41'' 


Quercetin 




2.60 ±0.04" 






Rutin 




10.42±0.33'' 







Each values represented in tables are means ± SD (N = 3}. 

Values with different letters (a, b, c, d, e, f, g, h) within same column were significantly different (P<0.05}. 
doi:l 0.1 371 /journal.pone.0096329.t004 
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Figure 2. Visual classification of HPLC chromatograms of Pyrola samples. 

dol:l 0.1 371 /joumal.pone.0096329.g002 



bath in the dark. Then, the absorbance of the samples was 
measured at 593 nm using the spectrophotometer. The trolox was 
used as the standard solution. The FRAP results were expressed in 
terms of micromoles trolox equivalent per g of dry weight (|tmol 
eq. trolox/g). All of the treatment groups were measured in 
triplicate. 

Data analysis 

The correlation coefficients of entire chromatographic patterns 
among samples were calculated, and the simulated mc-an 
chromatograms as well as characteristic peaks were generated 
using Computer Aided Similarity Evaluation (CASE) software as 
recommended by the Chinese Pharmacopoeia Committee. The 
software is used for evaluating similarities between different 
chromatograms based on the correlation coefficient (median) 
[17]. The results obtained in this study were calculated using the 
correlation coefficient unless otherwise specified. HCA and DA 
were performed using SPSS software (SPSS for Windows 18.0, 
SPSS Inc., USA) [18-20]. The 'average linkage between groups' 
method was applied and the cosine was selected as a measurement 
[21]. Three principal components obtained by principal compo- 
nent analysis were used to evaluate the similarities and differences 
among the tested samples [22]. 



Results and Discussion 

Validation of the HPLC procedure 

Methanol stock solutions containing tannin, hyperoside and 
quercetin were prepared and diluted to appropriate concentrations 
for the construction of calibration curves. Six concentrations of 
each analyte were injected in triplicate, and the calibration curves 
were constructed by plotting the peak areas under the curve versus 
the amount of the analytes. 

The precision of the analytical method was determined by 
assaying six replicates of the standard compounds (tannin, 
hyperoside and quercetin), and the relative standard deviations 
(RSDs) of the peak areas were estimated to be 0.38-1.16% (« = 6). 
The repeatability of the method was determined by injecting the 
same sample six times. The areas of the peaks were recorded, and 
the RSD of the areas varied from 1.13% to 2.49% (« = 6). To 
confirm the accuracy of the method, a recovery experiment was 
performed by mixing quantified samples with specific quantities of 
standard compounds. The average percentages of recovery of the 
three compounds ranged from 97.24±0.02% to 103.51 ±0.01%. 
In addition, the RSD varied from 1.13% to 2.20% {n-&). All the 
results demonstrated that the conditions of the analysis were 
repeatable and accurate (Table 2). 
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Figure 3. Dendrograms of hierarchical cluster analysis (HCA) 
for the tested samples of Pyrola. 

doi:l 0.1 371/journal.pone.0096329.g003 

Content of tannin, hyperoside and quercetin 

The contents of tannin, hyperoside and quercetin in three P)>rola 
species from different regions were analyzed by RP-HPLC. Based 
on a comparison of the retention times with those of the standards, 
peaks 1, 4 and 9 were identified as tannin, hyperoside and 
quercetin, respectively (Figure 1). The most abundant phenolic 
compounds in all samples were tannin and hyperoside, while 
quercetin was found at lower concentration (Table 3). The 
contents of tannin, hyperoside and quercetin varied within the 
ranges 9.77±0.01 to 34.75±0.05 mg/g DW, 0.34±0.01 to 
2.16±0.07 mg/g DW and ().062±0.004 to 0.147±0,007 mg/g 
DW, respectively. The content of tannin was significandy high in 
all tyrola samples, especially in P. calliantha. Hyperoside was the 
most abundant in P. renifolia. P. calliantha had higher contents of 
these three identified compounds among the three species, and P. 
calliantha from Left Banner of Alxa exhibited the highest content of 
tannin and quercetin among all samples. For each P. decorata 
sample from eight regions, the content of these three identified 
compounds varied greatly: the content of tannin, hyperoside and 
quercetin ranged from 9.77±0.01 to 16.48±0.08 mg/g DW, 
0.34±0.01 to 0.96±0.03 mg/g DW and 0.062±0.004 to 
0.080±0.006 mg/g DW, respectively. Tannin (16.48±0.08 mg/ 
g DW) and hyperoside (0.96 ±0.03 mg/g DW) were especially 
abundant in P. decorata from Fangyangsi, while quercetin 
(0.80±0.006 mg/g DW) was rich in P. decorata from Dianbing- 
chang. Recent experimental studies from our laboratory have 
shown that ecological factors, especially soil factors, contribute 
significandy to the content of tannin and quercetin in Pyrola 
samples, and this research also found that soil in Fangyangsi is 
richer in nutrients than other producing regions [8]. This could 
explain why P. decorata from Fangyangsi had high amounts of 
tannin and quercetin. 

Tannin, hypc-roside and quc-rcetin are all well-known human 
health antioxidants [23-25]. The present results showed that P. 
calliantha was rich in these active components, and could be used as 
a potential source of antioxidants for the food and drug industries. 
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Figure 4. Principal component analysis (PCA) of three principal 
components for Pyrola samples. 
doi:1 0.1 371/journal.pone.0096329.g004 

Total flavonoid content 

Flavonoids are the most common and widely distributed natural 
compounds, which are of great interest in nutrition and medicine 
Ixx'ausc of their antioxidant activity and possible pr()t(;('ti\-c effects 
as regards human health [26,27]. Comparing the three species, P. 
renifolia contained the highest average level of total flavonoid 
(27.26±0.69 mg/g DW), followed by P. decorata (22.54±6.61 mg/ 
g DW) and P. calliantha (18.09± 1.33 mg/g DW) (Table 3). As to P. 
decorata derived from eight regions, the highest and lowest total 
flavonoid contents were observed in P. decorata from Emei 
Mountain and Shanyang, respectively. The total flavonoid content 
of P. decorata from Emei Mountain was 37.82±0.54 mg/g DW, 
which was not only much higher than those of P. decorata samples 
from other regions, but also the highest of all samples. The total 
flavonoid content of P. decorata samples was in the following order: 
Emei Mountain > Fangyangsi > Xiabaiyun > Dianbingchang > 
Dadian > Ningshan > Lujuanliang > Shanyang. Based on these 
results, we may conclude that total flavonoid content varies gready 
among different species and among different geographical sources 
ot Pyrola samples. Zhang et al. reported that total flavonoid content 
of P. incarnata from different regions varied greatiy [1], which could 
support our findings reported above. 

DPPH radical scavenging activity 

DPPH assay has been widely used for the determination of 
antioxidant activity of pure antioxidant compounds as well as of 
different plant extracts, a lower DPPH ICso representing stronger 
antioxidant capacity [10]. For evaluation of antioxidant activity of 
tyrola, different samples were measured and their DPPH radical 
scavenging activity compared (Table 4). Among the three species, 
the highest radical scavenging activity was obtained for P. calliantha 
with the lowest average DPPH ICjo value of 9.66±2.09 |ig/ml, 
followed by P. renifolia (24.80±2.50 pg/ml) and P. decorata 
(37.11±8.16 |J,g/ml). DPPH radical scavenging activity of P. 
calliantha collected from Left Banner of Alxa (DPPH 
IC50 = 7.96±0.04 |tg/ml) was the highest among all samples, 
which was better than that of rutin. The DPPH radical scavenging 
activities of two other P. calliantha samples collected from 
Zhongshansi (DPPH IC50 = 9.02 ±0.03 Hg/ml) and Huzhu north- 
em Mountain (DPPH IC50 = 1 1-99±0.06 Hg/ml) were also high. 
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Table 7. Correlation coefficients between Individual chromatograms within a group and the group simulative mean 
chromatogram, and between the group simulative mean chromatograms. 
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It can be concluded that P. calliantha may be a valuable antioxidant 
natural resource. 

For P. decorata samples from eight regions, the DPPH ICr,,, 
ranged from 27.74±0.02 |Xg/ml to 50.33±0.01 ng/ml. The 
sample from Shanyang had the highest DPPH radical scavenging 
ability (DPPH IC50 = 27.74±0.02 |Xg/ml). These data indicated 
that the antioxidant activities of the same Pyrola species varied 
immensely from region to region. 

ABTS'^radical cation scavenging activity 

ABTS activity was quantified in terms of percentage inhibition 
of the ABTS** radical cation by antioxidants in each sample. The 
ABTS values of the thirteen samples were presented in Table 4. 
AH samples showed the capacity to neutralise the radical cation 
ABTS** and showed significant difference (P<0.05). Among the 
three species, the highest activity was obtained from P. calliantha 
with a value of 984.64 ±32. 48 |J,mol equiv. Trolox/g, followed by 
P. renifolia and P. decorata with values of 818. 96± 13.63 and 
696.78±67.94 |xmol equiv. Trolox/g, respectively. 

For P. decorata samples from eight regions, the ABTS values 
ranged from 612. 66± 1.26 to 798.23±7.33 |imol equiv. Trolox/g. 
The sample from Emei Mountain had the highest ABTS value 
(798.23 ±7.33 \imo\ equiv. Trolox/g). 

Ferric reducing power (FRAP) assay 

The FRAP assay evaluated the antioxidant properties of the 
samples based on their reducing ability. The values obtained from 
thirteen samples (Table 4) were significantly different (p<0.05). In 
details, P. calliantha possessed the highest antioxidant capacity with 
a FRAP value of 376.56±21.54 |j.mol equiv. Trolox/g, followed 
by P. renifolia and P. decorata with values of 325. 93 ±13.41 and 
268.48±42.08 |i,mol equiv. Trolox/g, respectively. 

As to P. decorata derived from eight regions, the FRAP value 
ranged from 219.64±4.66 to 31 1.45±2.04 |xmol equiv. Trolox/g. 
The sample from Shanyang had the highest FRAP value 
(31 1.45 ±2. 04 [xmol equiv. Trolox/g), followed by samples from 
Emei Mountain and Ningshan with values of 309. 12± 1.47 and 
306.45 ±3.94 |j,mol equiv. Trolox/g, respectively. 

The ABTS and FRAP values showed the same order of activity 
observed in the DPPH method. Based on these results, we 
concluded that Pyrola not only presented the highest free radical 
scavenge capacity but also the strongest reducing capacity. 
Meanwhile, the antioxidant capacity of Pyrola samples varied 
greatly of diHerent species of Pyrola and the same species from 
different regions of origin. 

In conclusion, through comprehensive analysis and evaluation, 
P. calliantha was selected as the most valuable species with the 
highest average content of total identified compounds (30.49 mg/ g 
DW). It also had the strongest average antioxidant capacity 
(DPPH IC5o = 9.66±2.09 |Xg/ml, ABTS va- 

lue = 984.64± 32.48 |Xmol equiv. Trolox/g and FRAP va- 
lue = 376.56±21.54 (imol equiv. Trolox/g), followed by P. renifolia 



and P. decorata. Emei Mountain was proved to be the most suitable 
region for producing P. decorata, because the highest content of 
combined phytochemicals (total identified compounds and total 
flavonoids) was observed in the Emei Mountain sample 
(50.85 mg/ g DW'). In addition, the antioxidant activity of Emei 
Mountain sample was the strongest among the eight P. decorata 
samples (DPPH IC50 = 29.91 ±0.03 ng/ml, ABTS va- 
lue = 798. 23± 7.33 ^mol equiv. Trolox/g and FRAP va- 
lue = 309. 12± 1.47 \ymo\ equiv. Trolox/g). 

HPLC fingerprint analysis 

Three Pyrola species from different regions were analyzed to 
develop a standard fingerprint under the established HPLC 
conditions. The simulated mean chromatogram was generated by 
CASE. Eleven common peaks were selected as characteristic 
peaks. Based on the comparison of the retention time and UV 
spectra with standard samples, peaks 1 , 4 and 9 were identified as 
tannin, hyperoside and quercetin, respectively (Figure 1). Peak 4 
(hyperoside) was chosen as the reference standard peak due to its 
peak area accounting for above 10% of the area of all peaks 
(Table 5). According to the analysis of chromatograms (Figure 1), 
we could roughly identify different species of Pjrola and the same 
species from different regions of origin. The CASE software was 
applied to evaluate the similarity of the chromatograms. The 
results showed that the correlation coefficients of similarity of 
chromatograms of thirteen samples ranged from 0.24 to 0.99 
(Table 6). The wide range of similarity revealed the differences of 
diverse samples. Furthermore, the similarity of the same species 
was higher, while that of different species was lower, these results 
also implying that the chromatograms were representative and 
associated with phytochemical constituents of Pyrola [28] . 

Hierarchical clustering analysis 

Prior to HCA, we compared the fingerprints visually and simply 
divided the samples into four distinct groups, namely A, B, C and 
D (Figure 2). Although it was possible to differentiate chromato- 
grams on the basis of visual comparison, this process was 
subjective and non-quantitative. HCA can provide a quantitative 
and objective analysis of fingerprints [20]. According to the 
chromatograms obtained from all samples, we selected eleven 
common peaks, using peak 4 as the reference standard peak, and 
then calculated their relative retention times and relative peak 
areas. A matrix was applied for HCA using SPSS software, which 
consisted of the number of the samples and the relative peak area 
of eleven constituents. A dendrogram was acquired using the 
average linkage between groups and the cosine method (Figure 3). 
At a clustering coefficient of 12.5, the thirteen samples can be 
divided into three groups. Groups A and D were merged into a 
new class, due to their similar chemical constituents, which 
revealed the subjectivity of visual classification. The three groups 
were renamed Gl, G2 and G3, which were consistent with the 
results of morphological classification. At a clustering coefficient of 
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Canonical Discriminant Functions 
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Figure 5. Canonical discrimination analysis (DA) of HPLC 
chromatograms for Pyrola samples. 
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5, the thirteen Pyrola samples were separated into six groups. In 
other words, HCA not only can distinguish different Pyrola 
species, but also can differentiate the same Pyrola species from 
different regions of origin. 

The correlation coefficients of each chromatogram for Gl, 
G2 and G3 corresponding to software-generated group simu- 
lated mean chromatograms are shown in Table 7. The 
correlation coefficients between these simulated mean chro- 
matograms are also shown in Table 7. The chromatograms 
within a particular group were generally consistent. Correlation 
coefficients for each chromatogram which were classified into a 
particular group to the software-generated group simulated 
mean chromatograms were higher than 0.90. However, the 
chromatograms within a particular group were markedly 
different from the chromatograms in different groups. Mean- 
while, the similarity of the software-generated group simulated 
mean chromatograms was lower than 0.75 and the difference 
between groups was evident [20] . Compared with similarities of 
the different groups, the ANOVN analysis revealed that the 
difference was significant between groups (P= 0.046). In 
principal component analysis plots (Figure 4), projected points 
for each group were localized in confined clusters, which can 
make a good explanation for the HCA results. 
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Discrimination analysis 

DA can be used to build a predictive model of the group 
membership based on observed characteristics in each case. This 
procedure can generate a discrimination function or a set of 
discrimination functions which can discriminate and classify the 
unknown membership on the basis of eigenvalues [20]. 

The relative peak areas of eleven common peaks were selected 
from the chromatograms. However, some peaks were meaning- 
less for establishing the discrimination function. Therefore, DA 
was applied to selected variables which were of value for the 
functions, and then the discrimination function was generated. 
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Two types of discrimination function were obtained using the 
SPSS software. 

Canonical discrimination function: 



Yi = 1.08x1-3.372x5+22.581x10 - 1.936 



Y2 = 0.804x1 + 2. 143x5 + 20.449xio - 3 .605 



Discrimination standard: 



7i>0 and Fi > - : Gi 



FkO and Fi < F2 : G2 



F2<0 and Fi < - F2 : G3 



Fisher's discrimination function: 



Gi= 7.912x1 -5.747x5 + 173.318x10 -17.897 



G2 = - 0.02x1 + 20.829x5 + 8.967xio - 16.572 



G3 = 1.899x1 - 1.210x5 + 36.232x10-1.951 

Discrimination standard: there were dependent variables, Gl, 
G2 and G,3, which denoted the samples from groups Gl, G2 and 
G3, respectively, x represented the independent variable (Table 8). 
We could obtain three functional values of each sample, and a 
sample was assigned to the group corresponding to the highest 
functional value. 

There were three variables which were selected to constitute the 
discrimination functions. When we needed to place an unknown 
sample, we inserted the relative peak areas of the three variables 
into the equations, and the unknown sample was grouped 
according to the discriminant standard value obtained. The well- 
resolved DA plots for the three groups are shown in Figure 5. 
Using the three most discriminating variables enabled tested 
samples belonging to groups Gl, G2 and G3 to be classified with 
100% accuracy. 

According to the HPLC fingerprint analysis, the chemical 
compositions vary gready among either different species or 
different geographical sources of tyrok samples. This variation 
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may lead to significant differences in effectiveness as functional 
foods and nutritional supplements. So it is important to evaluate 
the quality Pyrola. HPLC fingerprint analysis coupled with HCA 
and DA in the present research performed well for the quality 
evaluation of Pyrola samples. 

Conclusions 

The results from the present study indicated that all tested Pyrola 
samples significantly differed from each other in their phytochem- 
ical profiles and antioxidant property. Through comprehensive 
analysis and evaluation, P. calliantha was selected as the most 
valuable species, and the P. calliantha sample from Left banner of 
Alxa even possessed the strongest antioxidant activity among sdl 
the thirteen samples. Futhermore, Emei Mountain was proved to 
be the most suitable region for producing P. decorata. The results 
indicated that there were rich variations in Pyrola samples. 

Different chemical compositions may lead to significant 
differences in effectiveness as functional foods and nutritional 
supplements. In order to further evaluate the diversities and 
quality of Pyrola, HPLC fingerprint analysis coupled with HCA 
and DA were introduced to establish a simple, rapid and effective 
method for the accurate identification, classification and quality 
assessment of Pyrola. Thirteen samples were separated into three 
groups using HCA at a rescaled distance of 12.5, this separation 
being consistent with the results of morphological classification. 
Furthermore, the thirteen samples were separated into six groups 
at a rescaled distance of .5, which can be identified as the same 
species of Pyrola from different regions of origin. Two types of 
discriminant functions were generated using three selected 
predictor variables and the ratio of discrimination was 100%. In 
summary, this method can identify different species of Pyrola and 
the same species from different regions of origin. 

The present study provides meaningful information for the 
collection and application oi Pyrola in both healthcare and the food 
industry. Moreover, the results from this study could be used to 
compare and control the quality of Pyrola and other natural 
products prepared from them. It is well known that the species of 
Pyrola are very abundant; in this research we only studied three of 
them, so a large, systematic study of Pyrola from different sources is 
required in a future study. Additionally, it is understood that the 
DPPH scavenging capacity assay, ABTS*"^ radical cation scaveng- 
ing assay and ferric reducing power (FRAP) assay used in the 
present study may not fully reflect the antioxidative mechanisms of 
Pyrola in vivo and its health promotion properties in tht; human 
body, and further cellular and in vim studies of the biological 
activities are required. 
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